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Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field 
electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is 
demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow 
may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding 
the magnetic mean-field diffusivity in an electrically conducting turbulently moving compressible 
fluid. These phenomena occur if the irrotational part of the motion dominates the vortical part, 
the Peclet or magnetic Reynolds number is not too large, and, in addition, the variation of the 
flow pattern is slow. For both the passive scalar and the magnetic cases several further analytical 
results on mean-fleld diffusivities and related quantities found within the second-order correlation 
approximation are presented, as well as numerical results obtained by the test-fleld method, which 
applies independently of this approximation. Particular attention is paid to non-local and non- 
instantaneous connections between the turbulence-caused terms and the mean fields. Two examples 
of irrotational flows, in which interesting phenomena in the above sense occur, are investigated in 
detail. In particular, it is demonstrated that the decay of a mean scalar in a compressible fluid 
under the influence of these flows can be much slower than without any flow, and can be strongly 
influenced by the so-called memory effect, that is, the fact that the relevant mean-fleld coefficients 
depend on the decay rates themselves. 

PACS numbers: 52.65.Kj, 52.30.Cv, 52.35. Vd 



I. INTRODUCTION 

Many investigations of transport processes in turbu- 
lently moving fluids have been done in the framework of 
the mean- field concept. A simple example is the trans- 
port of a passive scalar quantity like the number den- 
sity of particles in a turbulent fluid [H-O- Another im- 
portant example is the magnetic-field transport in elec- 
trically conducting turbulent fluids. The widely elab- 
orated mean-field electrodynamics, or magnetofluiddy- 
namics, delivers in particular the basis of the mean-field 
theory of cosmic dynamos @, Q • 

The original equation governing the behavior of a pas- 
sive scalar in a fluid contains a diffusion term with a dif- 
fusion coefficient, say k. In the corresponding mean-field 
equation there appears, in the simple case of isotropic 
turbulence, the effective mean-field diffusivity k + in 
place of K, where Kt is determined by the turbulent mo- 
tion and therefore sometimes called "turbulent diffusiv- 
ity" . Likewise, the induction equation governing the 
magnetic field in an electrically conducting fluid contains 
a diffusion term with the magnetic diffusivity 77. In the 
mean-field induction equation there appears, again for 
isotropic turbulence, rj + rjt in place of 77, where rjt is 
again determined by the turbulent motion and, some- 
times, called the "turbulent magnetic diffusivity" . 

At first glance it seems plausible that turbulence en- 



hances the effective diffusion, corresponding to positive 
Kt and rit- In a compressible fluid, however, this is not 
always true. A counterexample for the magnetic case 
has been long known. Represent the velocity in the form 
It = V X i/j + V(/) by a vector potential ip satisfying 
the gauge condition V • i/j = and a scalar potential (f>. 
Let Uc, Ac, and Tc be a characteristic magnitude, length, 
and time, respectively, of the velocity field. Assume that 
the magnetic Reynolds number UcXcM is small compared 
to unity and that Tc considerably exceeds the free-decay 
time A^/?7 of a magnetic structure of size Ac. Then it 
turns out [7rtl0ll that 
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That is, negative rjt are well possible if the part of u 
determined by the potential 4> dominates. Then, the 
mean-fleld diffusivity is smaller than the molecular one. 
This surprising result deserves more thorough examina- 
tion, which is indeed one of the motivations behind this 
paper. Here it will be shown that a result analogous to 
dl]) applies to Kt, too. These results apply not only to tur- 
bulence in the narrow sense, but also to other kinds of 
random and even non-random (including steady) flows. 

Results of that kind might be of some interest for the 
turbulence in the interstellar medium. It is widely be- 
lieved that it is mostly driven by supernova explosions 
In this case the driving force, and so the flow, 
too, could have noticeable irrotational parts, i.e., parts 
that are described by gradients of potentials. However, 
when rotation or shear is important, or the Mach number 
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is close to or in excess of unity and the baroclinic effect 
present, vorticity production becomes progressively more 
important — even when the forcing of the flow is purely 



irrotational 15|. 



Another possible application of such results could be 
in studies of the very early Universe, where phase tran- 
sition bubbles are believed to be generated in connection 
with the electroweak phase transition [l^ ■ The rele- 
vant equation of state is that of an ultra-relativistic gas 
with constant sound speed c/\/3, where c is the speed of 
light. This is a barotropic equation of state, so the baro- 
clinic term vanishes. Hence, there is no obvious source of 
vorticity in the (non-relativistic) bulk motion inside these 
bubbles so that it should be essentially irrotational. This 
changes, however, if there is a magnetic field of signifi- 
cant strength, because the resulting Lorentz force is in 
general not a potential one. 

In Sec. II of this paper we give an outline of the mean- 
field theory of passive scalar transport, prove the passive- 
scalar version of relation ([Ij , and derive some further re- 
sults in the framework of the second-order correlation ap- 
proximation. We also give an analogous outline of mean- 
field electrodynamics and present some specific results. 
In Sec. Ill we formulate the mean-field concept for the 
case of non-local relationships between the turbulence- 
dependent terms in the mean-field equations and the 
mean fields, and we explain the test-field method for 
the determination of transport coefficients. In Sec. IV 
we then present analytical and numerical results for two 
simple models which refiect mean-field properties of ir- 
rotational fiows. Finally a discussion of our findings is 
given in Sec. V. 



II. OUTLINE OF MEAN-FIELD THEORIES 

A. Passive scalar transport 

Let us focus attention on passive scalars C which de- 
scribe the concentration of, e.g., dust or chemicals per 
unit volume of a fluid. We assume that C satisfies 



dtC + V ■ (UC) ~ V • (kVC) = , 



(2) 



where U is the fluid velocity and k a diffusion coeffi- 
cient, which in general depends on both the mass density 
g and the temperature T. In the case of an incompress- 
ible isothermal fiuid, Eq. ^ applies with k independent 
of position so that V • (kVC) turns into kAC. We want, 
however, to include compressible fluids, too. We may jus- 
tify ([2]), e.g., if C describes the concentration of an admix- 
ture of light particles in a compressible isothermal fiuid. 
The diffusion coefficient is then given by k = f{T)/g, 
with some function /; see [l^. Chap. §11, p. 39. We 
expect the validity of ([2]) with some dependency of k on 
g also in more general cases. In a flow of such a fluid, its 
density, even when uniform initially, will in general be- 
come position-dependent in the course of time. For the 



sake of simplicity we shall nevertheless ignore any con- 
sequence of inhomogeneous density. Hence, our results 
are applicable only for either the limited time interval or 
the limited velocity amplitude range for which the den- 
sity inhomogeneity is still negligible. Overcoming these 
limitations requires a theory which includes momentum 
and continuity equations and is beyond the scope of this 
paper. 

With this in mind we consider, in what follows, k al- 
ways as independent of position and replace ([2]) by 



dtC + V ■{UC)-kAC = 0. 



(3) 



We further assume that the fluid motion and there- 
fore also C show turbulent fluctuations, deflne mean 
quantities like C or by a proper averaging procedure 
which ensures the validity of the Reynolds rules, and put 
C = C + c and U — U + u. The evolution of C is then 
governed by 



dtC + V ■{UC) + g- kAC = 



with 



uc . 



For c we have 



dtc + V- {uC + Uc+ (uc)') - kAc = , 



(4) 



(5) 



(6) 



where (uc)' stands for uc — uc. Clearly, is a functional 
of u, U, and C in the sense that at a given point in 
space and time depends in general on it, f/ and C at 
other points, too. This functional is linear in C. 

Let us, for simplicity, consider the case U — and 
assume that u corresponds to homogeneous turbulence. 
Until further notice we adopt the assumption that C 
varies only weakly in space and time so that at a 
given point in space and time, can be simply represented 
as a function of C and its first spatial derivatives, i.e., 
VC, taken just at this point. We will refer to this as- 
sumption as "perfect scale separation". We may then 
conclude that 



dC 



(7) 



(c) 

with 7j and Kij being coefficients determined by m, 

which are independent of position.^ 

(c) — 
Clearly, 7^ gives the velocity of advection of C, and 

contribution to the total mean-field diffusivity 



^ Consequently, a spatially constant C is not influenced by u and 
thus stationary; nevertheless it causes a fluctuation c that is, 
for stationary u, again stationary. Hence, there is a non-trivial 
stationary solution of JSJ with constant average. For potential 
flows u = \/<j) it can be given explicitly as C = Co exp(0/K). We 
thank our referee for having made us aware of it. 
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tensor, which is then equal to K5ij + Hij. From ([7]) we 
conclude that 



^ = 7, 



dC 
dxi 



K 



(8) 



Of course, may be assumed to be symmetric in i and j. 

(c) 

For isotropic turbulence we have 7,^ = and — KtSij 
with some constant coefficient Kt so that 



and consequently 



-KtAC. 



(9) 



(10) 



Although we have defined u as the turbulent velocity, 
only its statistical symmetry properties like homogeneity 
or isotropy have in fact been utilized. Here and later in 
this paper, the term 'turbulence' should, accordingly, be 
understood in a wider sense, including random or even 
non-random flows with such properties. 

In specific calculations often the second-order correla- 
tion approximation (SOCA) is used. It consists in ne- 
glecting the term (uc)' in Eq. ^ for c, so that this 
equation turns into 



dtc - kAc = -V • (uC) . 



(11) 



The applicability of this approximation is restricted to 
not too large velocities u. If we characterize the veloc- 
ity field again by a typical magnitude and by typical 
length and time scales Ac and Tc , respectively, we may de- 
fine the parameter = A^/kTc, which gives the ratio of 
the free-decay time A^/k to t^, further, the Peclet num- 
ber Pe = UcXc/k and the Strouhal number St = Uc^c/Ac. 
Note that = Pe/St. A sufficient condition for the ap- 
plicability of SOCA in the case g^ ^ 1 reads St <C 1; in 
the case g^ <C 1 it reads Pe ^ 1. 



We further assume that u varies so slowly in time that 
we may consider it as independent of t. Finally we adopt 
SOCA so that (fTTj) applies. We may write it in the form 



A(kc- A) = 



(15) 



with 



X -X-G-2V<P-G + 4>G-x, A^ = (j). (16) 
From (fTSll we conclude that 



c — —X + Co . 

K 



(17) 



where cq is some constant. Calculating then ^ at x = 
we obtain 



1 



Ti = — (uiXk + 2u,d^/dxk)Gk ■ 
Due to the isotropy of the turbulence we have 



(18) 



UiXk = \u ■ x^ik , u.idf^ldxk = \u-V^bik- (19) 

Using (IT^ and and considering the homogeneity of 
the turbulence, we find 



Consequently we have 

^=-^iW-W)G. 
Comparing this with ([9]), we obtain 

«t = -^ii''^ - 0^) , 

6K 



(20) 



(21) 



(22) 



a result in full analogy to ([T]) . For an incompressible flow 
Kt can never be negative, while it can never be positive 
for an irrotational flow. 



B. Diffusivity in a special case 

Let us focus attention on homogeneous isotropic tur- 
bulence and determine Kt in a limiting case. Since Kt 
depends neither on C nor on position, we may choose 
simply C = G ■ X with a constant G so that VC = G, 
and consider at the end only a; = 0. We represent u in 
the form 



u = V xip + Vij), V-i>^Q, 



(12) 



by a vector potential ip and a scalar potential 0, and set 
further 



(13) 



where we utilized the freedom in defining the new vector 
potential x such that both ^p and (jj are now derived from 
this single quantity. We then have 



(14) 



C. Relations for transport coefficients 

We consider now homogeneous, but not necessarily 
isotropic turbulence and use a Fourier transformation of 
the form 

F{x, t)^ F{k, Lu) exp(ifc • x - iiot) d^k dio . (23) 



Further we adopt SOCA. Then standard derivations (see, 
e.g., 0) yield 



„(C) 



ifcfc 



Qik{k,uj) d^kduj 



Kk^ — ICJ 
Qij{k,uj) + Qji{k,uj) 



II ^ 2{Kk^-iuj) 

2K(Qik{k,Lo)kj + Qjk{k,uj)ki)kk\ , 

' d kdu! 



(24) 
(25) 



(k/c^ — iui)^ 
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where Qij{k,uj) is the Fourier transform of the correla- 
tion tensor Qiji^^r), defined by 



(26) 



Since Qij{^, r) is real, we have Qij{k, uj) — Q*j{~k, — cj), 
where the asterisk means complex conjugation. 

We recall here Bochner's theorem (see, e.c., 0, Chap. 
6), according to which, for any homogeneous turbulence, 
Qij{k,uj) is positive semi-definite, that is, 



Q,,{k,io)X,X* >0 



(27) 



for any complex vector X. 

Assume first incompressible turbulence, that is, 
V • It = 0. Then we have 



j kj — , Q ij ki — . 



(28) 



In this case, yields 7|*^'' — (even if the flow is not 
isotropic), and (f25|) turns into 



= ^ j j —g^—T^{Qij{k,Lu) + Qji{k,ijj)^ d^fcdw. 

(29) 

From (|27l) and we may conclude that k^- is positive 
semi-definite. If the flow is statistically isotropic we have 
i^ij = KtSij and we may conclude that is non-negative. 

Assume next irrotational turbulence, that is, u — 'Vcj) 
with any potential (j). Then we have 



Qijik^uj) = kikjR{k,uj) 



(30) 



with some real function R related to cj). Owing to (P7| . 
R must be non- negative (cf. 0, Chap. 6). With ([Ml) . 
(I25t and m we find 



7, 



(C) 



\kiR{k, uj)k^ 



d-^kduj 



II 



Kk^ — iuj 
(3Kfc^ -I- iuj)kikjR{k, uj) 
[nk"^ — iujY 



(31) 



d^fcdw. (32) 



For statistically isotropic flows R depends only via k on 
k. Hence, we obtain as expected 7, =0 and 



I 



(3kA:2 +iijj)k'^R{k,u:) „ 



(k/c^ — io;)^ 



d-'fcda;. 



(33) 



Assume now in addition that the variations of u in time 
are slow. Then R is markedly different from zero only for 
very small cj. Consequently, Kt is non-positive. 



D. Magnetic-field transport 

We now consider a magnetic field B in a homogeneous 
electrically conducting fluid and assume that it is gov- 
erned by 

dtB-V x{U xB)-r]V^B^Q, V • B = , (34) 



with U being again the velocity and rj the magnetic dif- 
fusivity of the fluid. Focusing attention on a turbulent 
situation, we define again mean fields, in particular B 
and U , and put B = B + b and U = U + u. Then we 
have 

dtB-'Vx{UxB + e)-rjV^B = 0, V -B = , (35) 
where 



£ = u X b 



(36) 



and 



dtb-V X {Uxb + uxB + {ux b)') - yyV^b = , 

(37) 



Here {u x b)' means uxb — ux b. The mean electromo- 
tive force £ due to the fluctuations u and b is a functional 
of u, U, and B, which is linear in B. 

Let us restrict ourselves again to U = 0. Assuming 
perfect scale separation, defined analogously to the pas- 
sive scalar case considered before, we may conclude that 



£i = aijBj - 7]ij{V X B)j - Cijk{VB] 



jk I 



(38) 



where (VB)^j, = ^{dBj/dxk + dBk/dxj). Here aij, r]ij 
and Cijk are quantities determined by u. (Instead of the 
traditional feyfe we use here rj^j = ^ bimnf^jmn and Cijk = 

^^{bijk + bikj)-Y 

In this context SOCA consists in dropping the term 
{uxb)' in ([SZ]) so that 



9tb-7/V^b = V X (m X B), V-b===0. 



(39) 



Sufficient conditions under which this applies are again 
analogous to those explained below (fTTj) . We have only 
to replace the parameter by = A^/yyrc and the 
Peclet number Pc by the magnetic Reynolds number 
Rm — UcAc/77. Note that = Rm/St. 

The relevant relations for aij, rjij, and Cijk, derived un- 
der SOCA for homogeneous turbulence, are given in the 
Appendix [X] Let us restrict ourselves here to homoge- 
neous non-helical turbulence. Then the correlation ten- 
sor Qij may not contain any pseudo-scalar or any other 
pseudo-quantity. As a consequence, the symmetric part 
of Qij and the antisymmetric part of 77^ are equal to zero, 
and we have 

(B) ff .^^ fe(fc,c.) + g,,(fc,^) ^3^^^ 

2 J J rjk/ — \u] 



^ In analogy to what was noted for the passive scalar case, sta- 
tionary solutions of l|34ll with constant mean parts of B are con- 
ceivable. 



5 



and 



1 f ff [25y 4; - {SikSji + SjkSii) 



2 JJ \ 2{r]k^ - iw) 

T] [2Sijkk - {hSjk + kjSjk)] h 



Qkiik,uj)d^kduj. (41) 



Moreover, Cijk is equal to zero. 

Consider now incompressible turbulence, for which 
(l28t applies. Then we have, even in the anisotropic case, 
7-'^^ — [see also Chap. 7.1, statement (i)]. Further- 
more, 



1 



[2SijSki ~ {SikSji + SjkSii)]Qkiik,uj) 3 



rjk'^ — iuj 



= - 1 1 — its — r. , 

(42) 

which, together with (j27l) . implies that rjij is positive 
semi-definite. If the turbulence is in addition isotropic, 
we have rjij = rjtSij with 



'7t 



1 /■ /■ Qkkjk.uj) 



d^fcdcj. 



(43) 



Like Kt, Tit also has to be non- negative [see also [tI, Chap. 
7.4, Eq. (7.47)]. 

Consider next irrotational turbulence, for which pop 
applies. Then, 



SB) 



Vij 



rjk'^ — lu) 

1 rr{T]k^ +iLj){k^Sij - kikj)R{k,uj) ^3 
2jj {r,k^ - iioY 



(44) 



d^fcdw. 



(45) 



If the variations of u in time are slow, R is markedly 
different from zero only for small lj. Then it can be 
readily shown that rjij is negative semi-definite. In the 

isotropic case, R depends, as already noted above, only 
via k on k. Therefore we have, independent of the time 
behavior of it, 7^ =0 and 

1 f f {r]k^ +iuj)R{k,Lu)k^ ^ 



Vt 



[rjk^ — iuj)^ 



d^fcdw. 



(46) 



If then the time-variations of u are slow, rjt has to be 
non-positive [see also 0, Chap. 7, Eq. (7.51)]. 



III. GENERALIZATIONS AND TEST-FIELD 
PROCEDURE 

A. Lack of scale separation 

In applications, the assumption of perfect scale sepa- 
ration, used so far, might be violated; see, e.g., [l^. We 
now relax it. Considering first again the passive scalar 
case we admit now a dependence of at a given point in 



space, on C (or its derivatives) at other points, that is, we 
admit a non-local connection between ^ and C. For the 
sake of simplicity, however, we assume until further no- 
tice that at a given time, is only connected with C (or 
its derivatives) at the same time, that is, we remain with 
an instantaneous connection between !F and C. Again, 
we restrict ourselves to homogeneous turbulence. 

We further assume here, again for simplicity, that 
mean fields are defined as averages over all x and y. 
Hence, they depend on z and t only. 

In what follows it is then sufficient to consider the z 
component of !F only. As a straightforward generaliza- 
tion of the relation for J^z contained in ([7]), with deriva- 
tives with respect to z only, we now write 



:F.{z,t)= / (ji^Hodz^ct) 



.(0 



dC{z-C,t) 
dz 



(47) 



with two functions 7!*^^ (C) and Kzz (C) i which are as- 
sumed to be symmetric in C, and with the integration 
being over all (In the case of inhomogeneous turbu- 
lence, 7i and Kzz would also depend on z.) With the 
specifications jz'^\C) = li'~^^S{C) and ^^^(C) = KzzS{C), 

(C) 

where 7i and Kzz on the right hand sides are under- 
stood as constants, we return just to the relation for jFz 
given by ([7]). Utilizing integrations by parts, we may 
rewrite (|T7l) as 



Tziz,t)= / r(C)C(z-C,t)dC (48) 



with 

r(C) = 7f' (C)-^^. (49) 

In what follows, it is convenient to work with a Fourier 
transformation defined by 

FiO^^ J F{k)eM^k(:)dk. (50) 

[Apart from the fact that here only a function of the sin- 
gle variable ( is considered, this definition differs from 
also by the factor l/27r on the right-hand side.] 
Equation (|48|) is then equivalent to 

^42, = 77- / i^{k)C{k, t) exp(ifcz) dfc , (51) 
2vr J 

and (I49p implies 

7f )(fc) ^ 5R(f(fc)) 
~^zz{k) = ~k-^^{t{k)). 



(52) 



For ^^^^ and Kzz on the right-hand sides of (O and ([5]) 
we have then 



7f)=7f'(0), nzz = k.M. 



(53) 
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Let us now admit that J^^, at a given time, depends 
on C (and its spatial derivatives) not only at this but 
also at earlier times. This non-instantaneous connection 
between T z and C can be described as a memory effect] 
see, e.g., [1^. We then have to generalize (|47| such that 



dC{z-C,t-T) 
dz 



dCdr (54) 



with 7z and symmetric in C, and equal to zero for 
T < 0; the integration is then over all C, and t > 0. It 
is straightforward to generalize the relations (|48l) to ((53)) 

in that sense. Then, Fourier transforms with respect to 

~ (c) 

C, and T occur, and 72 and kzz depend not only on k 
but also on an additional variable oj. 

The generalizations explained here can easily be ex- 
tended to the magnetic case discussed in Sec. Ill Dl Then, 
lz^\ Vxx, and rjyy occur as functions of ^, or of ^ and t, 
and their Fourier transforms jz , fjxx, and fjyy as func- 
tions of A: , or of A: and uj. 



1: 



B. Test-field procedure 



In Sec. |TT] we have presented results for quantities like 



or Kij which apply only under SOCA. As soon as we 

are able to solve equations like ©, e.g., numerically, we 

~ (c) 

may determine these quantities, or 7^ and kzz intro- 
duced in the preceding section, also beyond this approx- 
imation. A proper tool for that is the test-field method, 
first developed in mean-field electrodynamics [2l|, |2^ . 
We apply the ideas of this method here first to the passive 
scalar case. As in the preceding section we assume again 
that the mean fields are defined by averaging over all x 
and y and relax spatial scale separation, but ignore at 
first scale separation in time, that is, the memory effect. 

Suppose that we have solved ([6|) for two different test- 
fields C, say 



C = Co cos kz and C = Co sin kz 



(55) 



with given Co and fc, and calculated the corresponding 
J"z, say J-z{z) and J^zi^)- Specifying (ITH) to J^^ and 
and considering that, due to ([SH]) and the assumed 
symmetry of 7^ ^(C) and Kzz{C) in C: 



7f Ho cosfcCdC 



7f)(fc), 



we find 



KzziO cosfcCdC = kzz{k) , 



Co (t'z'^'' (fc) cos kz + Kzz (k) k sin kz) 
Co {ji^^ (k) sin kz — kzz {k) k cos kz) 



(56) 



(57) 



This in turn leads to 



1 ,^ 



^^^>{k) — — (J^zi^) cos kz + Fj^z) sinfcz) 
Co 

f^zz{k) — T^i-^zi^) sinfcz — J-liz) cos kz) . (58) 
Cok 

Note that, although constituents of the right-hand sides 
depend on z, the left-hand sides do not. 

If the memory effect is taken into account, Eq. ([6]) has 
to be solved with time-dependent test-fields C. Let us 
define such fields by multiplying the right-hand sides in 
([55]) by a factor e"^*. Integrate then the relevant equa- 
tions numerically with any initial condition until all tran- 
sient parts of the solutions have disappeared. For steady 
flows, the remaining solutions then show the same har- 
monic time variation as the test fields (approximately 
possible also for unsteady flows). That is, the same time- 
dependent factors occur on both sides of the equations 

analogous to (|57)) and can be removed. These equations 

(c) 

then allow the determination of 7^ (fc,w) and kzz{k,uj), 



that is, the Fourier transforms of 7; 

and kzz{k,LU 



(C) 



Of course, ^z'^\k,u} 



(C,t) and KzziCr). 
are in general com- 



plex. We may also replace the factor 



Aujt 



by 



with a 



as discussed in Sec. 
of quantities like 7] 



complex a. Instead of the Fourier transformation with 
respect to time we have then to use a Laplace transfor- 
mation. 

A test-fleld procedure, as described here for passive 

scalars, can also be established for the magnetic case 

II Dl It allows then the calculation 
jf — ' 

, Vxx, and rjyy or their Fourier or 
Laplace transforms. Such procedures have already been 
used elsewhere (e.g., [20l.[2a|). 

For the numerical computations presented below we 
use the Pencil Code [23], where the test-field meth- 
ods both for passive scalars and for magnetic fields have 
already been implemented All results presented in 
this paper have been obtained with a version of the code 
compatible with revision 16408. 



IV. EXAMPLES AND ILLUSTRATIONS 

A. Three-dimensional flow 

In an attempt to model properties of homogeneous 
isotropic irrotational turbulence we wish to consider first 
a steady potential flow. Thus, we choose 



(/> = 7^ cos ko{x + Xx) cos ko{y + Xy) cos fco(z 

Ko 



(59) 

Xz)- 

(60) 



Here, uq and ko are positive constants and Xx, Xy-i and 
Xz are understood as random phases. Of course, this 
steady flow must lead to growing inhomogeneities of the 
mass density. Therefore the applicability of our results 
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is restricted to a limited time range, see the discussion 
below ©. 

Starting from original fields C and B, which may de- 
pend on y, z, and t, and also on Xxi Xyt ^^'^ Xz, we 
define mean fields C and B by averaging over all x and 
y and, in addition, over Xz- Consequently, mean fields 
no longer depend on x, y, or Xz, but they may depend 
on z and t. Clearly, the Reynolds averaging rules apply 
exactly. For mean quantities determined by it only, aver- 
aging over X and y is equivalent to averaging over Xx a-nd 
Xy Therefore, such quantities can also be understood as 
averages over Xx, Xy, and Xz- Clearly, is independent 
of X, y, and also of z. 

From (l59l) and (l60l), we conclude that 




0.1 



1.0 10.0 
Pe or Rm 



1 /3 



and we define a wave number k{ of u by 
fcf — VSkf) . 



(61) 



(62) 



In what follows, nt and as well as kt and ijt, will be 
expressed in units of Kto and ijto, given by 



Kto = ■'yto 



3fcf 



(63) 



Furthermore, we define the Peclet number Pe and the 
magnetic Reynolds number Rm by 



T-t ^rms T-> ^rms 
Pe — ; — , Rm = 



nkf 



■qkf 



(64) 



Calculations in the framework of SOCA under the as- 
sumption of perfect scale separation yield 



= -KtoPe, 



(65) 



Clearly, Kt and ryt are non-positive. If scale separation is, 
in the sense of (I47L relaxed, we obtain 



FIG. 1: kt/hito and fjt/rito as functions of Pe and Rm, respec- 
tively, for the model given by ^ and k/kf = 1/10^/3 « 
0.06. The dotted line on the lower left gives the SOCA re- 
sult and the shaded area on the lower right marks the regime 
where the total diffusivities would become negative. 
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FIG. 2: Kt/^to versus k/k{ for some values of Pe. 



Kt(fc) 
^t(fc) 

fiv) 



-KtoPc f{k/k{), 
-?7to Rm/(fc/fcf) , 

1 + (2/3)1)2 



(66) 



In all following discussions we consider k as positive. Like 
Kt and Tjt, also Kt and fjt are negative as long as fc/fcf < 1. 

In what follows, we present results for the quantities 
Kt and fjt obtained by the test-field procedure described 
in Sec. IIIIBl utilizing numerical integrations of Eq. (|6]) 
for c or Eq. (|37|) for b. Averaging over Xz was performed 
by averaging over z. 

Figure[T]shows Kt / Kto as well as fjt /r]to for a small value 
of k/kf, at which these quantities should be very close to 
Kt/Kto and rjt/rjto, as functions of Pe and Rm, respec- 
tively. (These values could also be obtained with a test 
field that is independent of z and another one linear in 



z.) In agreement with the results presented in Sec.lTlland 
with (j65p . Kt and rjt are negative for not too large values 
of Pe and Rm, respectively. Remarkably the functions 
Kt(Pe) and 77t(Rni) coincide formally for small values of 
Pe and Rm only, but are otherwise clearly different from 
each other. In particular, % remains negative, at least 
for Rm < 70, while Kt becomes positive for Pe 2. The 
total diffusivities, 77 -I- r)t and k -I- Kti are always found to 
be positive. 

Figures [2] and [3] show examples of the dependence of 
kt/Kto and fjt/rjto on k/kf. Again, Kt and fjt with Pe = 
0.35 and Rm = 0.35, respectively, that is, in the validity 
range of SOCA, take coinciding negative values in the 
limit of small k/kf. However, Kt and fjt become positive 
for large values of k/kf, regardless of the values of Pe and 
Rm. The dependence of Kt on Pe and that of fjt on Rm 
are in general clearly different from each other. 
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FIG. 3: fjt/rjto versus k/kf for some values of Rm. 



The results regarding negative contributions of Kt to 
the mean-field diffusivity for passive scalars, or negative 
contributions of rjt to the magnetic mean-field diffusivity, 
have been found under the assumption that the velocity 
u is steady or varies only weakly in time. In order to 
see the influence of the variability of u we consider now 
a renovating flow. It is assumed that, during some time 
interval, a steady flow as given by (j60p with some values 
of Xxj Xy Xz exists, and likewise in the following 
interval, but with randomly changed Xx, Xy ^-nd Xz, and 
so forth. Hence, there is no correlation between the flows 
in the different intervals. It is further assumed that all 
intervals are equally long. Denoting their durations by 
r, we deflne now the dimensionless parameters 



(67) 



Steadiness of the velocity corresponds to = Qfj = 0. 

Figure |4] shows the dependency of Kt/z^toPe on and 
that of ryt/?7toR-in on for k/kf = l/\/3 ~ 0.6. We see 
that kt and fjt are no longer negative if and q^i exceed 
0.2 and 0.3, respectively. 



B. Plane wave-like flow 

With the idea of establishing a simple model reflect- 
ing features of homogeneous anisotropic turbulence, we 
remain with ((5^ . that is it = V(/), but replace (pO]) by 

4' ~ -j— cos {ko[sx -f z) — (jJot — x) ■ (68) 

If s = 0, the velocity u corresponds to a sound wave trav- 
eling in the z direction, with wavelength and frequency 
determined by feg and wq and with a phase angle x- We 
assume, for simplicity, fco > and wq > so that the 
wave travels in the positive z direction. If we admit non- 
zero values of s, the wave propagates no longer in the z 
direction, but in the direction of the vector (s, 0, 1). For 



FIG. 4: fi:t/«:toPe and i)t/r]toKra versus or g^, respectively, 
for a renovating flow with UrmsfctT ~ 5.3 and fc/fcf = l/\/3 « 
0.6. In the calculations, r was held constant. Consequently, 
= Pe/5.3 and = Rm/5.3. 



loq — Q the velocity u does not depend on time, that is, 
we have a "frozen-in" wave. 

Similarly to the preceding example we deflne mean 
fields here by averaging over all x and y and, if the orig- 
inal field depends on also over x- Then, mean fields 
may depend only on z and t. Again, the Reynolds averag- 
ing rules apply exactly. If an original field is determined 
by u only and s is unequal to zero, averaging over x is 
equivalent to averaging over x- 

Instead of (ICTI) we now have 



1 



and instead of and ([55)) we put 

k{ ^ ko, Kto = Vto = 



(69) 



(70) 



and we define Pc and Rm again according to ([M)) . Finally 
we set 

Due to the definition of mean fields, which implies that 
C and B do not depend on x and y, we also have (V x 
B)z = 0. In addition, we assume that = 0. 

Adopting SOCA and assuming again perfect scale sep- 
aration, we find 



7; 



and 



(C) 

'(s,g) 



f^zz — 

"UXX = Tjyy 



UimsPe g{s,q^) , 
Ui.msRni5(s,q^) , 

(1 + s^)q 
{l + s^y+q^ ' 



-KtoPe h{s,qK.) , 
-T]to Rm h{s, qr,) , 



(72) 



(73) 
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h{s,q) 



(l+s2)((l + s2)2_3q2) 



1.00 



((1 



5^)2 + 



In the case = qr; = 0, that is, for frozen-in waves, 

the 7i and jz vanish. This is due to the fact that, 
then, there is no preference for the positive or negative 

(C) 

z direction. For 7^ 0, however, 7^ is positive so that 
C is advected in the positive z direction. Further, Kzz 
is negative for not too large ^k, but it becomes positive 
for larger q^ . This applies analogously with g,, , 7i , and 
^xx — Vyv' 

Relaxing perfect spatial scale separation and assuming 
that C and B vary in time as exp(crt) with a real a, we 
find further 

7(C') =u„„3Pe g{s,qf^,k/k[,a/Kki) 

if'' ^ Uy.^s'Rm g{s,qr,,k/k[,a/T]kf) (74) 

g / 1 + s2 + „ 

g[s, q,v,w) = - "2 

^ V {{I + v)^ + + w) +g2 

1 + s2 - w \ 



((1 - t;)2 + s2 + w) + g2 



and 



fizz = - Kto Pe h{s, q^, k/k{, aj Kkf) 

fjxx = Vyy = - ritoRmh{s,qr,,k/kf,a/r]kf) (75) 



1 / (l + s2+^;)((l+^;)2+s2_^^^ 

2y \ {(l + vy + s^ + w) +q-^ 

_ ^ (1 + - -^)((1 - Vf + 5^ + \ 

((1 - i;)2 + s2 + w)^ + g2 j 

The factor (5 is in general equal unity but equal to zero 
ifl — w = s = w = g = 0, that is, if the following 
denominator vanishes. All coefficients i'f"\ if\ ^^^\ ^^f\ 
Kij, rjij, kij, and 77^ , which are not explicitly mentioned, 
are equal to zero. 

Numerical calculations of 7!'^'' and iizz as well as 
and fjxx = fjyy by the test-field method, without restric- 
tion to SOCA, have been carried out with k/kf ~ 0.1 
and some specific values of Pe and Rm. Only cases with 
s were included, for which the x x averages 
are equivalent. Hence the standard test-field procedure 
with horizontal averages could be employed. Figure [5] 
shows these quantities for s = 0.01 as functions of g^ or 
q,j. The results for Pc = 0.1 and Rm = 0.1 are in good 
agreement with our SOCA calculations, that is, ([71]) and 
([75)) . Those for higher Pe and Rm clearly deviate from 
them. Interestingly, the dependence of jz and kzz on 
Pe is always the same as those of jz and fj^x or fjyy on 
Rm. 
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FIG. 5: Dependence of 7i'^V'^rmsPe and kzz/ntoPe on 

as well as that of 7i^'/urmsRni and fj^x/vto^T^ on Qn for 
the model given by Eq. ^ for k/kf = 0.1, s = 0.01, 
CT = and three values of Pe or Rm, respectively. Solid 
lines give the SOCA results, symbols the values obtained by 
the test-field method; dots correspond to 7z'^V'^rmsPe and 
Kzz/wtoPe, open circles to ^i^^ /umns'Rm, and fj^x /ritoRrn. 
Clearly, ji'~^' /uunsPe and Kzz/KtoPe coincide completely with 
'yi^'/urmsRm and ?7a;a;/j?toRin if Pe and Rm coincide. 



A remarkable feature of the SOCA results ([75)1 for kzz, 
and also for fj^x = Vyy, is that these quantities show sin- 
gularities at /c/fcf = 1 if = (7,, = s = (7 = 0. Neverthe- 
less they are well defined at this point; kzz/n — Pe^/4 
and fixx/v = Vyy/V = Rm^/4 at k/kf = 1. 

Let us, in what follows, focus attention on passive 
scalars only. Consider a mean scalar of the form C = 
F{t) cos kz with F being positive. For q^ = its time be- 
havior is exclusively determined by the quantity k + kzz- 
Clearly C is bound to decay if k -I- kzz > 0. Now con- 
sider the dependence of kzz on k/kf for g^ = s = a = 
depicted in Fig. [S] If fc/fcf is smaller than but close to 
unity, kzz may, even for small Pe, take arbitrarily large 
negative values, and n + kzz becomes negative. This will 
then lead to a growth of the modulus of C. Of course, 
this conclusion is drawn from a result obtained under the 
SOCA and may hence be questionable. Indeed, the suf- 
ficient condition for the applicability of SOCA given so 
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/C„+K<0 
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FIG. 6: Dependence of kzz/f^ with Pe = 0.707 and = 
a = on k/kf. Solid lines as well as the cross at k/k{ — 1 
result from analytic SOCA calculations with s = 0. Symbols 
give numerical results obtained with Eq. © using s = 0.01; 
filled circles: full equation, open circles: SOCA, term (nc)' 
dropped. Note the second "resonance" at fc/fcf = 2. 



far, Pe ^ 1, has been derived for k/kf <C 1 only. If, by 
contrast, fc/fcf « 1, we find, when comparing the terms 
V • (uc)' and kAc in ([SJ under the assumption that c 
is dominated by contributions with wavenumbers fc + fcf 
and k — k{, for — a = and s <^ 1 the more stringent 
condition 



Pe< 3 



(1- fc/fcf)^ + 
1 - fc/fcf + s2 



(76) 



It supports the doubts in the above conclusion concerning 
the growth of the modulus of C. 

The aforementioned SOCA calculations for s = have 
been extended by the inclusion of fourth-order terms in 
It, that is, in Pe. Apart from some quantitative changes 
of Kzz/k in the neighborhood of fc/fcf = 1, which occur 
with larger Pe, a new singularity emerges at fc/fcf — 2. 
As can be seen in Fig. [51 the numerical (non-SOCA) cal- 
culations with s — 0.01 reflect this feature, too. They 
also give indications of a further resonance at fc/fcf — 3 
(not shown). 

Figure [7] shows k^z for steady test fields (that is, a — 
0) with fc/fcf = 0.9 and s = 0.01 as a function of Pe. 
The results clearly deviate from those obtained by the 
SOCA as soon as Pe exceeds, say, 0.2. Considering that 
K + > is equivalent to —kzz/Kto < 1/Pe, Figure [7] 
tells us further that k + kzz becomes very small with 
growing Pe, but suggests that it remains positive. We 
may suppose that the modulus of the considered C never 
grows but its decay becomes very slow for large Pe. For 
example, for Pe = 1 we expect that A = — (k + kzz)k'^ ~ 
— O.lftfc^, that is, the decay of C should be about ten 
times slower than in the absence of any motion. 

In these considerations, however, the memory effect, 
that is, the dependence of the value of kzz, relevant for 



FIG. 7: Dependence of kzz/n%o on Pe for k/ki — 0.9 and 
s = 0.01 cr = g„ = 0. The dashed hue gives the SOCA 
result and the shaded area marks the range where the total 
diffusivity would become negative. 
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FIG. 8: Dependence A(cr) = -[{k + k^^ik, a)]k^ for Pe = 1.0, 
k/kf = 9/10, s = 0.01, and = 0. The curve representing 
\{a)/Kk'^ intersects the dash-dotted line A = o" at X{a)/Kk^ = 
—0.005, which is shown more clearly in the inset. 



the decay of C, on the decay rate A itself, has been ig- 
nored. As explained in Sec. IIII B[ we have to include 
this dependence by using time-dependent test fields. Let 
us assume that they vary as exp(at) but consider a first 
as independent of A. Then kzz and A ~ — (k + kzz)k'^ 
occur as functions of a. Figure [SI obtained by test-field 
calculations, shows this dependence of A on cr. If we 
then identify a with A we find, as indicated in Fig. [51 
A w -0.005Kfc2. That is, the decay of C is about 200 
times slower than in the absence of any motion. 

In order to check this surprising result, we perform two- 
dimensional direct numerical simulations based on Eq. 
^ with a flow given by using fco = 10 fci, where fci = 
2tt / Lz is the smallest wave number in the z direction with 
extent Lz- Our computational domain is periodic in both 
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FIG. 9: Direct simulation showing the decay of the ampli- 
tude F of C{z) in units of its initial value Fq. Parameters like 
in Fig. |8] The inset shows the time dependence of the growth 
rate A, leveling off at A ~ —0.005 K,k^ after some initial ad- 
justment time. 
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FIG. 10: (Color online) Snapshot of C{x,z) for the simula- 
tion shown in Figure |9] at tKk^ — 17. Bright (yellow) shades 
indicate positive values and dark (blue) shades negative val- 
ues. The lower panel shows C{x, z) with the correct aspect 
ratio of the box. 



directions. We choose = 10 Lz so as to accommodate 
the variation in the x direction with wavenumber fco* 
and s = 0.01. The initial condition is C = Cocosfcz, 
with Co > and k = 9ki. We use 128^ mesh points and 
choose Pe = 1, which is clearly beyond the applicability 
of SOCA; cf. Fig. H 

As we expect that C = F{t) cosfcz, we have identified 
the maximum of the a;-average of C with respect to z 



with F and determined the growth rate by calculating 
first its instantaneous value X{t) — dlnF/d<; see Fig. [HI 
It turns out that the average of A over the time interval, 
in which it is approximately constant is in excellent agree- 
ment with the test-field result A ~ — 0.005Kfc^, described 
above. The snapshot in Fig. [TUl shows that C{x, z) varies 
mainly in the z direction with the dominant wave num- 
ber k = 9fcf/10. In units of nk"^ , the free-decay rate of a 
mode with this wave number would be 0.01, or 0.012, if 
the variation in x is taken into account. Note, however, 
that the dominant constituent of C belongs, by virtue of 
its X dependence, to the fluctuating field c and that the 
actual decay rate of the mean field C is at least two times 
smaller than the given free-decay rate. The fluctuations 
are not decaying freely, but follow the mean field, and 
hence adopt its decay rate. In this particular case, the 
rms values of the fluctuations exceed those of the mean 
field by a factor of 14. 

The question could be raised as to whether a resonance 
effect in the above sense can also occur for solenoidal 
flows. Numerical experiments with the (stationary) ABC 
flow (for its definition see, e.g., [l^l) indicate clearly that 
the decay of C is always accelerated in the presence of 
this flow, irrespective of the value of k/kf. 



V. CONCLUSIONS 

In this work we have shown that the turbulent diffu- 
sivity Kt for the concentration of a passive scalar in a po- 
tential flow can be negative at low Peclet numbers. This 
result is analogous to an earlier finding for the turbulent 
magnetic diffusivity ryt in such a flow at low magnetic 
Reynolds numbers, originally derived in the context of 
astrophysical dynamo theory. The numerical calculations 
presented in this paper confirm Eq. ([T]) quantitatively for 
an irrotational flow. We have not yet considered the case 
of the combined action of solenoidal and irrotational flows 
where the question arises of how strong the solenoidal 
part has to be to render the turbulent diffusivities posi- 
tive. Our calculations also show that negative values of Kt 
do not occur for larger Peclet numbers, whereas negative 
77t may well exist for moderate to large Reynolds num- 
bers. In neither case have negative turbulent diffusivities 
yet been seen in laboratory experiments. Nevertheless, 
for possible physical applications of our results one may 
think of microfluidic devices 25], in which the flow can 
be compressible [1^ and the Peclet number small. 

In addition to the condition of small Peclet and mag- 
netic Reynolds numbers, there are also the requirements 
of good scale separation and of slow temporal variations 
of the flow. If these requirements are not obeyed, Kt and 
rjt are no longer necessarily negative - even at small val- 
ues of Peclet and magnetic Reynolds numbers. This may 
be the reason why a reduction of the effective diffusivity 
has never been seen in physically meaningful compress- 
ible flows and why Eq. ([1]) is virtually unknown in the 
turbulence community. In fact, previous attempts to ver- 
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ify this equation in simulations have failed because of the 
fact that the time dependence has been too vigorous in 
those flows prj . 

The spatial structure of the flow does not appear to 
be critical for obtaining a reduction of the effective diffu- 
sivity. Even in a nearly one-dimensional flow, turbulent 
diffusivities can become negative. However, in that case 
there are two new effects. First, if the underlying flow 
pattern displays propagating wave motions, there can be 
transport of the mean scalar in the direction of wave 
propagation - even in the absence of any mean mate- 
rial motion. Again, this effect may have applications to 
microfluidic devices. Second, the wavenumber character- 
istics display a singular behavior under SOCA, but even 
beyond SOCA there can be a dramatic slow-down of the 
decay by factors of several hundreds compared with the 
molecular values. This result is completely unexpected 
because no such behavior has ever been seen for any other 
turbulent transport process. Furthermore, the memory 
effect proves to be markedly important in such cases, 
so the common assumption of an instantaneous relation 
between the mean flux of the scalar and its mean con- 
centration or the mean electromotive force and the mean 
magnetic field breaks down near the singularity. 

In addition to finding out more about possible appli- 
cations of the turbulent transport phenomena discussed 
above, it would be natural to study the possibility of 
similar processes for the turbulent transport of other 
quantities including momentum and heat or other active 
scalars. Further, a complementary effort to determine 
the transport coefficients for turbulent irrotational flows 
numerically would be of high interest, the more as there 
are no simple analytical results available. Supernova- 
driven turbulence in the interstellar medium would of 
course be the most suggestive application. 

Clearly, both analytical and numerical approaches us- 
ing the test-field method proved to be invaluable in that 
they are able to predict unexpected phenomena that can 
then also be verified using direct numerical simulations 
and in future hopefully also laboratory experiments. 



Appendix A: Relations for aij, riij, and Cijk 

Under SOCA we may derive, for homogeneous turbu- 
lence, 



K'^iimkj - eujkm) ^'"l^^'^^ d^kduj , (Al) 
rik^ — luj 



, , , r r X s '^Vi^ijklkjyi kikmSjl) 



rjk'^ — iu) 



Qlm{k,Uj) 3 

X — -r^ — r-^ d^k doj , 



(A2) 



-2r]- 



'^^ijnnkjkf^ (^f^ijnjkk ^imkkj')ki 



rjk^ — icj 

, Qmn j3 

rjk-^ — iuj 



d^kduj. 



(A3) 
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